EPJ manuscript No. 

(will be inserted by the editor) 



The effects of grain shape and frustration in a granular column 
near jamming 



o 

(N 

C 
3 



J.M. Luck 1 and A. Mehta 2 

1 Institut de Physique Theorique, IPhT, CEA Saclay, and URA 2306, CNRS, 91191 Gif-sur-Yvette cedex, France, 
e-mail: jean-marc.luck@cea.fr 

2 S.N. Bose National Centre for Basic Sciences, Block JD, Sector 3, Salt Lake, Calcutta 700098, India, 
e-mail: anita@bose .res . in 



I 

-a 
c 

o 
o 



> 

m 
m 

o 
o 



X 



Abstract. We investigate the full phase diagram of a column of grains near jamming, as a function of 
varying levels of frustration. Frustration is modelled by the effect of two opposing fields on a grain, due 
respectively to grains above and below it. The resulting four dynamical regimes (ballistic, logarithmic, 
activated and glassy) are characterised by means of the jamming time of zero-temperature dynamics, and 
of the statistics of attractors reached by the latter. Shape effects are most pronounced in the cases of strong 
and weak frustration, and essentially disappear around a mean-field point. 

PACS. 64. 60. My Metastable phases - 45.70. Cc Static sandpiles; granular compaction - 45.70. Vn Granular 
models of complex systems; traffic flow - 64.70.Q- Theory and modeling of the glass transition 



1 Introduction 

One of the reasons why heterogeneities are intrinsic to 
granular media is the absence of thermal motion; spatial 
and temporal behaviour has no reason to equilibrate, so 
that structures as well as time tracks which are far out 
of equilibrium can remain juxtaposed in the same system. 
Clearly, this leads to spatially and dynamically hetero- 
geneous behaviour. It is, however, only recently that re- 
search efforts in this context have focused on heterogeneity 
(see [1] for a recent review). Examples of static spatial het- 
erogeneity in granular systems include bridges |2|3|4j and 
force chains [5] . More generally, dynamical heterogeneities 
have attracted a lot of attention recently, both in granu- 
lar matter (see [5]) and in other systems such as glasses 
or colloids (see [?])• 

Spatiotemporal heterogeneity takes place when differ- 
ent parts of a system have diverse dynamical behaviour 
characteristic of their location. The first indications of 
such heterogeneity in a vibrated granular medium were 
found in the experiments of Reference [5] , the findings of 
which indicated that both the average density, as well as 
density fluctuations, varied strongly throughout a shaken 
box of grains as a function of depth. Computer simula- 
tion and theoretical results [5] reproduced this behaviour, 
predicting additionally that the mean density was an in- 
creasing function of depth, and that density fluctuations 
were largest in the middle of the box for time windows rel- 
evant to experiment. These results additionally suggested 
that the phase behaviour within the box was very hetero- 
geneous; ballistic behaviour was expected near the top, 
activated in the middle and glassy behaviour at its base. 



The theoretical model [1011 1112113] on the basis of 
which the latter predictions were made is as follows: gra- 
ins in a column are able to orient themselves in one of 
two possible ways, corresponding to 'ordered' and 'disor- 
dered'. When a grain is in its disordered orientation, space 
is wasted: a void of size e is created, which characterises 
the shape of the grain. Rational and irrational values of e 
correspond to regular and irregular grain shapes respec- 
tively. Despite the simplicity of this 'aspect ratio' formu- 
lation of shape effects, recent work [14] has shown that 
it may be used to characterise a rich variety of granular 
shapes. 

The presence of gravity is included in the model by a 
depth-dependent local frequency, such that lower (more 
weight-bearing) grains move more slowly than the less 
burdened upper grains. In the jamming limit, voids are 
minimized; accordingly, the definition of a ground state in 
the model is one that locally minimizes the voids ratio [5] . 
However, since the orientation that minimizes voids with 
respect to grains above a given grain is not typically that 
which fulfils the same function for those below itself, this 
naturally generates frustration. In the model, this is repre- 
sented by the effect of two oppositely directed fields whose 
relative strengths are modulated by a coupling constant g. 



In earlier work [10111112113) , the focus was on the effect 
of shape; the behaviour of the model for typical rational 
and irrational values of the shape parameter e was ex- 
plored for g = [11112] and then in the g — > limit [L3l. 
Here we complete the analysis by looking at the effect of 
varying the coupling constant g. This is equivalent to vary- 



2 



J.M. Luck, A. Mehta: The effects of grain shape and frustration in a granular column near jamming 



ing the frustration, and as will be shown, has wide-ranging 
effects on the behaviour of the model. 

The plan of the paper is as follows. The definition of 
the model is recalled in Section 2. Section 3 is devoted to 
the statics of the model, i.e., the number and the nature of 
its ground states. We then address the properties of zero- 
temperature dynamics. Section 4 contains an analysis of 
the dynamical phase diagram and of the behaviour of the 
jamming time in the various regimes, whereas the statis- 
tics of attractors is investigated in Section 5. In Section 6 
we conclude with a brief discussion of our findings. 

2 The model 

In its most complete form, the model |13j consists of a 
finite column of N grains, labeled by their depth n = 
1, . . . , N. Each grain has an orientation variable o~ n = ±1. 
Grain n is referred to as up or ordered when o~ n = +1 
and down or disordered when er n = —1. Disordered orien- 
tations generate voids and waste space, whereas ordered 
ones do not. Implicit in this description is the effect of 
shape, which is most easily understood in terms of the 
rectangular grains of aspect ratio a considered in [10 . 
Grains aligned along their long edges (length 1) result 
in a fully packed column, whereas those perched on their 
short edges (length a < 1) leave voids of size e = 1 — a. The 
horizontal orientation is thus ordered, and the vertical one 
disordered. Such a two-state model is clearly an approx- 
imation; we lump the effects of all possible void spaces 
created by disordered orientations of arbitrarily shaped 
grains into one disordered (vertical) orientation, and make 
a similar approximation for the ordered (horizontal) ori- 
entation. 

The N binary variables {a n = ±1} define the 2 N 
configurations of the system. We consider the following 
continuous-time stochastic dynamics which do not obey 
detailed balance. Grain orientations are updated with the 
Markovian rates 

f w(a n = +1 a n = -1) = e -( A "+ ff ")^, . , 
\ w(a n = -1 a n = +1) = e-<^- ff »)/^ [ZA) 

where 

• r is a dimensionless vibration intensity, referred to as 
temperature. 

• A„ is the activation energy of grain n, which we take to 
be proportional to its depth: 

A„ = (2.2) 

sdyn 

The dynamical length £dyn is the depth beyond which gra- 
ins are frozen out by the sheer weight of grains above 
them. Thus, the frequency of response of a grain n falls 
off exponentially with its depth: 

0„ = c- A "/ r = e-"/^". (2.3) 

• H n is the local ordering field felt by grain n, which is 
determined by all the other grains, both above and be- 
low n. The effect of the upper grains is assumed to be 



uniform. The back-propagation from grains below a given 
grain cannot, of course, be similarly uniform. We assume 
for simplicity that upward constraints are exponentially 
damped, with a characteristic length £i nt , the interaction 
length. We thus write 

H n = h n +gj n , (2.4) 

where the uniform effect h n of grains above n (m = 1, 
. . ., n — 1) and the non- uniform effect j n of grains below n 
(m = n + 1, . . . , N) are given by 

n-l N 

K=J2f( a rn), 3n= £ /(*») e^""^ , 
m,—l m— n+1 

(2.5) 

whereas g is a positive coupling constant. 

Furthermore, both components h n and j n of the total 
local field H n acting on grain n depend on every grain ori- 
entation <7 m = ±1 through the same function, the 'shape 
factor' /(cr m ), where 

/(a) = £zI_£±± ff= J* ^ ° "J. (2.6) 
J y ' 2 2 [-1 if (7 = 11. v ' 

The parameter e can be thought of as the size of a 
typical void space for a grain of a particular shape, with 
rational and irrational values of e corresponding to regular 
and irregular grain shapes respectively. There is an exact 
symmetry between models with e and 1/e, so that the 
shape parameter can be restricted to < e < 1. 

Putting all of this together, we see that that the con- 
tribution of an ordered grain to the local field is (negative) 
unity, while that of a disordered grain is a void space of 
magnitude e; the latter is clearly a function of granular 
shape, hence the name given to the shape factor f(a m ). 
The minus sign in the ordered case ensures that the con- 
tribution of an ordered grain to the excess void space H n 
is less than that of a disordered grain, as it ought to be; 
more importantly, this says that a void space is destroyed 
every time a grain aligns in an ordered fashion relative to 
its neighbours. 

The dynamics of the model involves the often conflict- 
ing contributions of the opposing local fields h„ and j n . 
In turn, these comprise all the terms f(o~ m ) which take 
values e or — 1 according to (|2.6|) . for all the other gra- 
ins m in the column. As mentioned above, this represents 
a simple-minded way of incorporating frustration into the 
model. 

In the following, we use the notation 

*int -e- 1 ^'". (2.7) 
We mention the following recursion relations: 

fl n = h n -l + /(<7„_x), j n = Eint (f((7 n +l) + jn+l) > 

with hi = Jn = 0, which provide a fast algorithm to 
evaluate the local fields. 

To sum up, the model parameters are the number 
of grains N, the shape parameter e, the coupling con- 
stant g, and the interaction and dynamical lengths £; n t 
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and £dyn- Previous work has been devoted to investiga- 
tions of zero-temperature static (number and structure of 
ground states) and dynamic (recovery of ground states as 
attractors) properties of the model in several special cases 
of interest: the directed model (g = 0) for £dyn = 00 [H]> 
the directed model for general £d yn [H], and the weak- 
coupling regime (g<C 1) for e = 1 |13) . 

In this work we aim at giving an overall picture of 
zero-temperature properties of the model all over its pa- 
rameter space, with an emphasis on their dependence on 
the coupling constant g and on the lengths £j n t and £dyn- 
As some of the features of the model are different for ra- 
tional and irrational values of the shape parameter e, we 
shall use e — 1 as our prototypical rational number, and 
the (small) golden mean 



f 



'5-1 



0.618033 



(2.9) 



<P 2 

as our prototypical irrational number. 

3 Statics: ground states 

The rules (|2.1I) simplify as follows in the zero-temperature 
limit (r -+ 0): 



w(a n = -1 ->■ cr n = +1) 
w(<Jn = +1 -> cr„ = -1) 



= e 2H ^ r -> 



00 if H n > 0, 
if H n < 0. 

(3.1) 

It is therefore natural to define a ground state as a con- 
figuration where the orientation of every grain is aligned 
along its local field |l()lllll2ll3j : 



sign H„ 



+ 1 if H n > 0, 
-1 if H n < 0. 



(3.2) 



We start with two special cases where the analysis of 
ground states is simpler. 



3.1 Directed model (g = 0) 

The ground states of the directed model have been inves- 
tigated in [11112] . In that case, the local field H n = h n 
acting on grain n only depends on the grains above n. 
The expression (|3.2j) therefore boils down to the following 
recursion relation: 



h n > => cr„ = +1, h n+ i = h n — 1, 
h n < =+ er„ = -1, hn+i = h n + £, 



with initial values h\ = 0, and o\ = +1 for definiteness. 
In a ground state, all the local fields h n lie in the range 



1 <h n <e. 



(3.4) 



The boundedness of the local fields implies that all the 
ground states have the same mean orientation (a), such 
that (f(o-)) = 0, hence 



g- 1 

7TT' 



(3.5) 



up to fluctuations which become negligible for large sys- 
tems. 

The number and the nature of ground states depend 
on whether e is rational or irrational. 

If the shape parameter e is irrational, the recursion 
formula (|3.3p implies that all the local fields h n are non- 
zero (except hi = 0). A unique quasiperiodic ground state 
is thus generated. Had we made the initial choice o~\ = 
— 1, we would have obtained the same configuration, up 
to a permutation of the two uppermost grains, so that 
the model has in all two ground states. For the golden 
mean (|2.9[) . the ground-state grain configuration is given 
by a Fibonacci sequence: 



-+- 



If the shape parameter e is rational, i.e., 

P 



(3.6) 



(3.7) 



in irreducible form (p and q mutual primes), some of the 
local fields h n generated by the recursion (|3.3p vanish. 
The corresponding grain orientations a n remain unspeci- 
fied. This orientational indeterminacy occurs at points of 
perfect packing, such that n — 1 is a multiple of the pe- 
riod p + q. The model therefore has extensively degenerate 
ground states. Every one of them is a random sequence of 
two well-defined patterns of length p + q, such that each 
pattern contains p up and q down grains. Defining the 
static (configurational) entropy per grain as 



S = 



InM 



A' 



N 



(3.8) 



where A/jv is the number of ground states of a system 
consisting of N grains, we have therefore 



E = 



In 2 

p + q 



(3.9) 



in the limit of a large system. The simplest of all rational 
values is e = 1, i.e., p = q = 1. In this symmetric case, 
we have f(o~) = —a, so that both orientations play sym- 
metric roles. The ground states are all the dimerised con- 
figurations, made of the patterns H — and — K The static 
entropy here assumes its maximal value S = (ln2)/2. 



(3.3) 3.2 Mean-field point {g = 1, £ int = 00) 



This situation is the complete opposite of the previous one. 
Upper and lower grains have equal weights, so that both 
components h n and j„ of the local field add up to give 

H n = r]-f(a n ), (3.10) 

where we have introduced the mean field 

N 

V=J2 ^( <7 ™) = eN ~ -N+=eN-(e + 1)N+, (3.11) 

771 — 1 
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with N + and N~ — N — N + being respectively the num- 
bers of up grains and of down grains in the configuration. 
The above expression shows that 77 is a global measure of 
excess void space [5] in the system. This globality results 
from the exact cancellation of fluctuations in the local void 
space corresponding to the competing fields h n and j n ; 
mean-field behaviour thus replaces the local fluctuations 
of the general case. In this limit, the model resembles the 
one studied in [TU], one of the earliest building blocks of 
the present model. We shall comment on further analogies 
between both models in due course. 
The condition Q3.2p thus reads 



sign(?7 - /On))- 



(3.12) 



It is fulfilled for all grains n as soon as the mean field lies 
in the range — 1 < rj < e. These inequalities amount to 
saying that N + takes a well-defined ground-state value: 



Int 



sN + 1 
e + l 



(3.13) 



where Int(x), the integer part of x, is the largest integer 
less than or equal to x. The result (|3.5[) is recovered in the 
limit of a large system. 

At the mean-field point, the ground states are all the 
configurations consisting of iV+ up grains and N~ = N — 
N^ s down grains. The number of ground states is therefore 



AO 



N 



N 
N+ 



so that the static entropy per grain reads 
£ = ln(e + 1) — Ins. 

^ ' £ + 1 

This result will be illustrated in Figure [5] 



(3.14) 



(3.15) 



3.3 General case 

We now turn to the general case. The behaviour of the 
model turns out to be dictated mainly by the coupling 
constant g, whereas the effect of the other parameters e 
or £i n t is less pronounced. The difference between rational 
and irrational values of e, which is the most salient feature 
of the directed model, manifests itself most in the weak- 
coupling and strong-coupling regimes (j < 1 and g ^ 
1). The overall picture, already sketched in [13], is the 
following. 

If the shape parameter e is irrational, the static en- 
tropy rises continuously from zero and behaves linearly at 
weak coupling: 



£^Ag ( 3 «1) 



(3.16) 



where the amplitude A depends on e and £;„t. The mean- 
ing of this result [13] is that generic ground states consist 
of quasiperiodic patches whose typical length diverges as 
h(g) ~ 1/g at weak coupling. The entropy then smoothly 



increases as a function of g, reaches a maximum around 
the mean-field coupling g = 1, and smoothly falls off to 
zero for j> 1. 

If the shape parameter e is rational, the static entropy 
stays equal to its value (|3.9[) over a whole range < g < g s , 
where g s (with 's' for static) is the static threshold. As 
already underlined in [13] . evaluating g s is a non-trivial 
task in general. The simplest situation is for s = 1 and 
an even number of grains (N — 2K). The ground states 
in the directed case {g = 0), and by continuity at weak 
enough coupling, are the 2 K dimerised ones. It can then be 
argued, thinking along the lines of [13] , that the first non- 
dimerised ground states which appear upon increasing g 

are + H ( — \-) K ~ 2 and h + (H — ) K ~ 2 , and that the 

relevant grain orientation for their stability is the second 
one (0-2)- Consider the first configuration for concreteness. 



We have ho 



-1 and j 2 = x int 



■T" X: r 



h int ' 



X; 



~ v "" ' ' The above configurations become ground 



int ^int 

states for H 2 
thus obtain 



fh 



h 2 +932 > 0, i.e., g > g s = \h 2 \/j 2 - We 



1 + Xint 



X int (l + 2x int 

and for an infinitely large system 



ix int x int / 



,9s 



1 + Zint 



a; int (l + 2x int + 2a; int ) 



(3.17) 



(3.18) 



Figure [T] shows a plot of the static entropy £ against g 
for a column of = 20 grains with £j nt = 10, for both 
shape parameters e = 1 and e = 1/<P. The exact number of 
ground states is obtained by means of a full enumeration 
of the 2 20 configurations. The distinction between rational 
and irrational e is visible at weak coupling. For s = 1 
(our prototype of a rational), S remains equal to £ = 
(ln2)/2 0.346573 in the whole range < g < g s , where 
the static threshold reads g s w 0.491651, including the 
finite-size effect of (|3.1T|) . For e — 1/<P (our prototype of 
an irrational) , £ rises continuously from the minimal value 
(In2)/A^, i.e., essentially zero, up to a finite-size effect due 
to the existence of two ground states. 

Besides this, the static entropy has a weak depen- 
dence on the shape parameter e. Both for rational and 
irrational e, the entropy has a smooth maximum around 
the mean-field point (g = 1), and it falls off smoothly at 
large g. The increase in the static entropy as the coupling g 
increases from to 1 reflects the progressive decrease of 
the amount of order in the ground states. For an irra- 
tional e, this corresponds to a shortening of the coherence 
length h(g) with increasing g. Finally, the weak depen- 
dence of the entropy on e near its maximum, i.e., near the 
mean- field point, is an early indication that shape depen- 
dence is increasingly lost as the model become more and 
more mean-ficld-likc. 



A digression on the effects of shape 

As mentioned above, increasing the value of g corresponds 
to increasing the frustration. As soon as the coupling con- 
stant g takes appreciable values, the evolution of ordering 
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1 „ 2 3 



Fig. 1. Plot of the static entropy £ against the coupling 
constant g for a column of N = 20 grains with £i n t = 10 for 
both shape parameters e = 1 and e — 1/<1>. For e = 1, the first 
symbol is at the threshold g s « 0.491651 (see dSTTTJ)). 



due to compaction no longer proceeds in a top-down fash- 
ion, as it did in the directed (g = 0) model (10111112] . 
We outline here what we might expect for the dynamical 
behaviour in the case of general g, before making specific 
calculations. 

Looking at zero-temperature dynamics allows us to 
get a flavour of the ordering behaviour. Do grains re- 
trieve their ground states in the limit of zero perturba- 
tion, and if so, how? In the case of g = 0, the effect of the 
shape parameter e is maximal on the zero-temperature dy- 
namics [10111112] : irregularly shaped grains with irrational 
e order ballistically fast into their unique quasiperiodic 
ground state, while regularly shaped grains never retrieve 
any of their many ground states, manifesting density fluc- 
tuations instead, due to the presence of many sites where 
h n = [11112] . Similar behaviour has been observed in 
experiments [8] on vibrated granular beds, from which we 
can infer that the effects of frustration, if any, are negli- 
gible compared to the intensity of vibration in the exper- 
iment under consideration. 

For a low but non-zero coupling g, the 'reverse' field j n 
induces some frustration as it begins to impose its own or- 
der in an upward direction, competing with the downward 
ordering due to h n . In the case of irrational e, one might 
imagine that the effect of the reverse ordering would break 
up the unique quasiperiodic ground state corresponding 
to g = 0; while in the case of rational e, even the tiniest 
amount of frustration ensures that points of zero field are 
not constantly generated and re-generated, and hence that 
the density fluctuations of the 3 = case disappear. The 
only possible logical culmination of ordering from both the 
top and the bottom of the column might be to 'fix' the 
points of zero field; in the case of e — 1 this corresponds 
to a dimerisation of the ground states. These were indeed 
the findings of earlier work [13] . 

Increasing frustration beyond these values might in a 
general sense lead to the breakdown of even this level 
of order. While the details would depend on parameters 
like £i nt and £dyn, we mention a few likely outcomes. For 
both weak (g <C 1) and extremely strong (g ^> 1) cou- 



plings, we might expect the prevalence of very similar 
order, although propagating in opposite directions (top- 
down in the first, and upwards in the second case). At 
the mean- field point (g — 1, £i n t = £dyn = oo), where 
the fluctuations in both ordering fields cancel each other, 
the system of grains is totally uncorrected, leading to 
a situation similar to that explored in [10]. The gradual 
replacement of the ordering by individual grains by the 
ordering of granular clusters might be expected to occur 
for g in the vicinity of the mean- field point (g = 1); shape 
effects are therefore expected to be minimal here. 

The results of the next two sections will bear out some 
of these speculations. 

4 Dynamics: jamming time 

The rules for zero-temperature dynamics are defined as 
follows [11112113] , The uppermost grain is kept fixed to 

<Ji=+l. (4.1) 

The other grains (n — 2, . . . , N) are selected at a rate 
given by (|2.3p . Once a grain is selected, its orientation a n 
is aligned along the local field H n according to the deter- 
ministic rule 

cr„^signiJ„, (4.2) 

provided the local field H n does not vanish. The choice 
of boundary condition (|4.ip is motivated by the fact that 
the strongest component in the local field H n is typically 
the long-ranged component h n which propagates via grav- 
ity. As a consequence, interactions propagate downwards 
in general, so that it is natural to impose a boundary 
condition at the top of the column. We assume that the 
column is prepared in a random state, where each grain is 
oriented at random (o~ n — ±1 with equal probabilities for 
allra>2). 

The above rule is well-defined for a non-zero coupling 
constant <?, because the local fields H n do not vanish in 
general. The zero-temperature dynamics thus defined le- 
ads to metastability. A finite column of N grains is even- 
tually driven to an absorbing configuration or attractor, 
in a finite jamming time T. This attractor is necessarily 
one of the ground states described earlier, i.e., a configu- 
ration where every orientation <r„ is aligned with H n . Let 
us emphasize that imposing a restrictive boundary condi- 
tion, i.e., fixing one of the spins (see (|4.1jl ). is necessary to 
have metastability in the above sense. Without such re- 
striction, zero-temperature dynamics would drive a finite 
system to a fluctuating steady state. 

In the present context, 'metastable state', 'attractor' 
and 'ground state' are therefore essentially synonymous. 
Arbitrary initial conditions can lead to any one of the 
ground states being reached. They are however fragile, 
in the sense that a slightly different initial condition or 
stochastic history leads to another attractor being reached 
in general. This fragility [T3] of metastable states is one of 
the characteristics of granular media jTB] . 

Along the lines of [13] , we will focus on two aspects of 
zero-temperature dynamics, namely the statistics of the 
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jamming time (in this Section) and that of the attractors 
(in Section 5). The jamming time is doubly random, as 
it depends both on the initial configuration of the system 
and on its whole stochastic history. Just as for statics, we 
begin with a few special cases. 



4.1 Directed model (g = 0) 

The dynamical behaviour of the directed column has been 
studied at length |11I12| . It depends qualitatively on whe- 
ther e is rational or irrational. 

If the shape parameter e is irrational, its unique qua- 
siperiodic ground state is reached by ballistic coarsening. 
An upper layer of the column is ordered, whose thickness 
grows linearly with time (the ballistic phase will be de- 
scribed more thoroughly in Section 4.3). The dependence 
of the corresponding velocity V on e has been investigated 
in fT2]. 

If the shape parameter e — p/q is rational, the local 
field h n may vanish whenever n — 1 is a multiple of the 
period p + q. It is therefore natural to complete the dy- 
namical rule P~2"j) as [11112] : 



± with prob. 1/2 



if 
if 
if 



K > o, 

h n = 0, 
h n < 0. 



(4.3) 



These dynamical rules do not drive the system to any 
of its ground states. There are always grains whose local 
fields h n vanish. The column reaches a non-trivial fluc- 
tuating steady state, investigated in |12) . which exhibits 
anomalous roughening: the fluctuations in the local field 
grow with a power law, as (h%) ~ ti 2//3 . 



4.2 Mean-field point (g = 1, £mt = oo) 

It has been shown in Section 3.2 that the statics of the 
model is of a mean-field type when g = 1 and £i n t = oo. 

This property extends to the dynamics in the £dyn = 
oo limit, where activation energies are negligible, so that 
grains are sampled uniformly and the effect of gravity is 
lost. In this limit, the dynamical rule (|4.2I) indeed becomes 



sign(r?- /(ovO), 



(4.4) 



where rj is the mean field introduced in (|3.11[) . Thus every 
grain orientation a n is updated to +1 (resp. —1), with 
unit rate, irrespective of its position n, as long as rj > e 
(resp. rj < —1). This rule can be recast as the following 
effective dynamics for the number N + of up spins: 

N+ <N+ =>■ N+ -> N+ + 1 with rate N - N+, 
N+ >N+ N+ -> N+ - 1 with rate N+ . 

(4.5) 

The dynamics stop as soon as N + reaches the value N+ 
(see (|3.13|) ). i.e., when the system reaches a ground state. 

Mean-field zero-temperature dynamics are fast, in the 
strong sense that the jamming time is microscopic. More 



precisely, consider s < 1 for definiteness, so that iV+ < 
N/2. For a random initial configuration, characterised by 
N + « N/2, the mean jamming time can be shown to read 
(see e.g. [17] ) 



(T) 



N/2 

E 

N+=N+ 



1 

N^ 



In. 



2s 



(4.6) 



for a large system. The jamming time is indeed found to 
be microscopic. This makes good physical sense, since the 
system is fully uncorrelated, and all the grains are simul- 
taneously mobile. 

Furthermore, as N + is the only non-trivial dynamical 
variable, it is clear that all the ground states are reached 
with uniform probability. In other words, anticipating the 
discussion of Section 5, the mean-field dynamics of our 
column model (g = 1, £j nt = £dyn = oo) is one of the 
rare instances where Edwards' flatness hypothesis [T5] can 
be shown to be exactly valid. A similar result has been 
established in the context of the ageing dynamics of mean- 
field spin-glass models [19]. 



Role of a finite £d yn 

We now study the model at its static mean-field point, but 
with generic dynamics defined by a finite value of the dy- 
namical length £dyn- This situation is of interest because 
it is both simple (the updating rule is still given by (|4.4p ). 
and non-trivial (grains are not selected uniformly any- 
more). Grain n is indeed updated at a rate given by (|2.3[) . 
so that upper grains are more mobile than lower ones. 
This case is physically similar to that of a column of non- 
interacting grains in the presence of gravity |10| . 

Consider again e < 1 for concreteness. For a random 
initial configuration, the typical number of grains to be 
flipped from + to — reads 



N f 



N 
~2 



N, 



2(1 + e) 



■N. 



(4.7) 



It is worth considering first the slow regime (N S> 
fdyn)- In this situation, the grains to be flipped are essen- 
tially the Nf uppermost + grains, which occupy an upper 
layer of depth 2Nf in a random initial state. As a conse- 
quence, the jamming time scales as T ~ exp(2ATf/£dyn)) 



i.e., 



T ~ exp 



1-e N 

1 + £ Cdyn 



(4.8) 



The jamming time therefore grows exponentially with N. 
The prefactor is proportional to that entering the asymp- 
totic orientation (13.51) . As a consequence, (|4.8[) does not 
apply to the symmetric situation (e = 1), where the jam- 
ming time is not exponentially large in N in this regime. 
Large jamming times testify that the retrieval of ground 
states is not easy; they are likely to correspond to lower 
dynamical entropies. The e-dependent prefactor suggests 
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that shape effects would be strongly related to the ab- 
sence of Edwards flatness, as will be verified below. Fur- 
thermore, the typical orientation profile of an attractor 
takes the form of a discontinuous step: 



-1 (n<2JV f ), 

(2N{ <n<N). 



(4.9) 



We define the zero-temperature dynamical entropy per 
grain as 

s = in q(c), ( 4 - 10 ) 

c 

where the sum runs over the attractors C, and where Q(C) 
is the probability that the dynamics drive the system into 
attractor number C, starting from a random initial con- 
dition. The above picture of jamming in the slow regime 
leads to the estimate NS « (N - 2jV f ) In 2, as the differ- 
ence N — 2N{ is an estimate of the number of the lower 
grains which do not move during the history of the col- 
umn. The dynamical entropy therefore reads 



5' 



2 s 



1 



In 2. 



(4.11) 



This quantity is smaller than the static entropy S, given 
by (|3.15l) . Both entropies indeed only coincide at the ex- 
tremal values e = (where £ = S = 0) and e = 1 (where 
S = S = In 2). The entropy difference is maximal for 
e = 1/4, where it equals S - S = ln(5/4) « 0.223143. 
Figure [5] presents a comparison between both entropies. 
This slow mean-field regime is one of the rare cases where 
the violation of Edwards' flatness can be turned into a 
quantitative estimate. Another zero-temperature exam- 
ple is provided by kinetically constrained one-dimensional 
spin models [20] . 




Fig. 2. Comparison between static and dynamical zero- 
temperature entropies against the shape parameter e in the 
range < e < 1, in the slow regime of the mean-field model 
(g — 1, £mt — oo, N ^> Cdyn)- Upper full curve: static en- 
tropy E (see (|3.15[) '). Lower dashed curve: dynamical entropy S 
(see (EH)). 



For generic values of the ratio iV/£dyn, the attractor 
statistics vary continuously between the uniform case of 



mean- field dynamics (for N <C £dyn) and the non- uniform 
case of the slow regime, described above (for N 3> £dyn)- 
This continuous dependence is best visualised by the ori- 
entation profile (a n ) of the attractors. Figure [3] shows nu- 
merical data for the orientation profile with e = \/<P and 
N = 200. Each dataset is obtained by averaging over 10 6 
different stochastic histories with different initial configu- 
rations. The data demonstrate a continuous crossover be- 
tween a uniform profile at the mean value (|3.5p . i.e., (a) = 
-l/<£ 3 sa -0.236067 (for N < £ dyn ) and the discontinu- 
ous step profile (|4.9[) of the slow regime (for N ^> £dyn)- 
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Fig. 3. Plot of the orientation profile (a„) of the attractors 
against depth n at the static mean- field point (g = 1, £j n t = oo) 
for e = 1/& and N — 200. Symbols: data for several values 
of Cdyn- Dashed lines: limiting uniform and step profiles, re- 
spectively corresponding to JV < £dyn and N ^> £dyn- 



4.3 Full dynamical phase diagram 

We now turn to the zero-temperature dynamics of our 
model for generic parameter values. We expect that the 
phase diagram of the model can be roughly divided into 
three regions: 

(a) the weak-coupling regime (g <C 1), already explored 
in earlier work for e = 1 |13j : a strong uniform field h n is 
the dominant effect, leading to a correspondingly strong 
dependence on the shape parameter e, 

(b) the neighbourhood of the mean-field point (g = 1), 
where ordering proceeds with fewer local constraints, so 
that shape dependence is increasingly lost, 

(c) the strong-coupling end (j> 1), where a strong frus- 
trating field j n is the dominant effect; this might be ex- 
pected to lead to the return of a strong dependence on the 
shape parameter e. 

For the time being, we restrict our study to the case 
where £ d yn = oo, so that grains are sampled uniformly 
by the dynamics. The dynamical properties of the model 
are mainly dictated by the coupling constant g and the 
interaction length £j nt , with a less pronounced dependence 
on the shape parameter e, except in the weak-coupling 
regime. The main features of the model are summarised 
in the dynamical phase diagram shown in Figure 21 This 
picture will be made more precise in the following. 
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Fig. 4. Sketch of the phase diagram of the model in the g— 
£int plane for a rational ((R), left) and irrational ((I), right) 
shape parameter e. Symbol: mean-field point. Curves: critical 
lines fj n t = Cint,c(3;e)- B: ballistic phase. A: weak-coupling 
activated phase. A': strong-coupling activated phase. 



At this point, it is worth emphasising that the dy- 
namical phase diagram of the model depends on the pre- 
cise definition of zero-temperature dynamics, including 
the boundary condition (|4.1[) . In particular, had we cho- 
sen to fix the bottommost spin (ov/v) instead of the upper- 
most one (cti), we would have obtained somewhat different 
phase boundaries; this would be true especially for large g, 
with a finite limiting £i n t, c along the g — oo axis mirroring 
that which exists currently for g = 0. Wc will return to 
this point in Section 6. 

The scenario observed for generic values of g qualita- 
tively follows that of the weak-coupling (g -C 1) regime 
for e = 1, investigated in [T3J. The model is in a ballis- 
tic phase if the interaction length is small (£j nt < £int,c)j 
and in an activated phase if the interaction length is large 
(£int > 6nt,c)- The critical value £- m t,c of the interaction 
length £; nt depends strongly on g and weakly on e. Let us 
start by reviewing the main characteristics of both phases 
and of the crossover between them, which have been ana- 
lyzed in [15] . 

• Ballistic phase (£i n t < £i a t,c)- In this phase, zero-tempe- 
rature dynamics propagate order into the system from the 
top down [11112113] . More precisely, if we define the thick- 
ness L(t) of the upper ordered layer of the column as the 
depth of the uppermost grain which is not aligned with its 
local field, the ballistic phase is characterised by a linear 
growth of the mean thickness: 



(L(t)) » Vt, 



(4.12) 



where V is the ballistic velocity. Fluctuations around this 
mean behaviour are due to diffusion. As a consequence, 
for a finite system of N grains, the jamming time grows 
linearly with N: 

CO « y, (4.13) 

up to relatively negligible fluctuations, so that the reduced 
variance of the jamming time, 



K T = 



varT _ (T 2 ) 



(T) (T)' 



-1, 



(4.14) 



is of order 1/N. 

• Activated phase (£; n t > £mt,c)- In this phase, zero-tempe- 
rature dynamics do not proceed in any ordered way. The 



system explores its configuration space more or less uni- 
formly, until it meets one of its ground states by chance. 
This picture is that of an activated phenomenon. An expo- 
nential growth of the mean jamming time with the column 
size results: 

(T) ~ e aN , (4.15) 

at least for very large N, where a is the effective reduced 
activation energy per grain. In other words, the system 
has to cross an extensive entropic barrier, whose height 
grows asymptotically as aN, in order to reach a ground 
state. This also suggests an exponential distribution of 
jamming times, so that the reduced variance asymptotes 
to Kt = 1. 

• Diffusive crossover (£i nt « £int,cj- The crossover between 
ballistic and activated behaviour has been shown [T3] to be 
described by a simple effective model. The thickness L(f) 
of the ordered layer was treated as a collective coordi- 
nate, and its dynamics modelled by biased Brownian mo- 
tion (see Appendix B of [T3J). In the ballistic regime, the 
downward propagation of the layer is helped by the dom- 
inant effect of the field h n , while in the activated regime, 
it is hindered by the dominant effect of the field j n . The 
crossover thus corresponds to the point £i n t = £i n t,c where 
the effects of the two fields are neutralised. The behaviour 
right at the crossover is dictated by the presence of a 
diffusive critical point. Observables obey finite-size scal- 
ing laws involving the scaling variable z = aX + ft, with 
X = Ar(£j nt — £int,c)> whereas a and (3 are non-universal 
numbers. For instance, the mean jamming time (T) and 
its reduced variance Kt obey 

(T)k^F(z), K T ^G(z), (4.16) 

where D is an effective diffusion constant, whereas F and 
G are known universal scaling functions. The reduced vari- 
ance of the jamming time takes the non-trivial universal 
value Kt — G(0) — 2/3 right at the critical point, i.e., 
for z = or £i nt ~ £i»t,c — (3/{aN). In practice the 1/N 
correction is negligible for N/£- m t lC > 10. Measuring Kt 
thus provides an efficient way of exploring the dynamical 
phase diagram. 

We now turn to the presentation and discussion of ac- 
tual numerical data. Figure [5] shows a plot of the ballistic 
velocity V against the coupling constant g, for both shape 
parameters e = 1 and e = 1/<P, and for two values of 
the interaction length (£i nt = 3 and 5) deep in the ballis- 
tic phase, i.e., much smaller than £i n t, c - For the irrational 
shape parameter e = 1/<P, the velocity V departs contin- 
uously from its value in the directed model (g — 0), i.e., 
V w 2.58 [12]. It increases steadily and reaches its max- 
imal value at or near the mean-field coupling [g = 1). If 
the shape parameter e is rational, the velocity assumes a 
constant value over a whole range < g < gd, where g^ 
(with 'd' for dynamical) is the dynamical threshold. Along 
the lines of [13], the latter can be shown to read 



9d 



1 ~ Xint 



1 



a; in t max(p, q) max(p, q) 



(4.17) 
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for an arbitrary rational value e = p/q of the shape pa- 
rameter. For £ = 1, this result reads 



1 



3d 



Xint = e i/fo 



1. 



•^int 



(4.18) 



This dynamical threshold is always smaller than its static 
counterpart (see (|3.18p ). as it should be. Both thresholds 
will be shown in Figure [5] 




0.2 0.4 0.6 0.8 1.0 



Fig. 5. Plot of the ballistic velocity V against the coupling 
constant g for several values of the parameters £i n t and e. 



to fall off as £int^ 2 (because of the law of large numbers). 
The formula (I4.6[) shows that the corresponding times also 



1/2 



to 



fall off as £ int . This brings in an extra factor of £ 
the velocity, providing a plausible mechanism for the ex- 
ponent 3/2 in P~T9]1 , 
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Fig. 6. Plot of the ballistic velocity V at the mean-field cou- 
pling (g = 1), divided by the interaction length £i nt , against 
£ ; Y t 2 , for both shape parameters e = 1 and e = 1/$. The dashed 
line has slope 2.2. 



The increase of V with g to a maximum near the 
mean-field point, observed for both irrational and ratio- 
nal e, seems apparently paradoxical, as increasing g after 
all increases the effects of frustration. The dynamics at the 
mean-field static point (g = 1) are however expected to be 
rather fast, irrespective of shape, so that the velocity can 
be expected to be both large and independent of e near 
the mean-field point. We will return to this paradox later, 
when we investigate the scaling properties of the jamming 
time near the mean-field point. 

Right at the mean-field coupling (g = 1), the ballis- 
tic velocity V is observed to depend strongly on £i nt and 
weakly on e, along the lines of the above discussion. It 
diverges faster than linearly at large £i n t. In order to in- 
vestigate this divergence, we present in Figure |5] a plot 
of the ratio V/£i n t against dot' for both shape parame- 
ters e = 1 and e = 1/<P. The largest velocities reached are 
V ~ 1000 for £j nt = 50; velocities larger than this are hard 
to measure accurately. The observed linear behaviour of 
both datasets with equal slopes suggest the scaling be- 
haviour 

V(g = l)^A^, (4.19) 

the prefactor A « 2.2 being seemingly independent of 
shape parameter. The exponent 3/2 of the growth of the 
velocity, as mean-field dynamics are approached, can be 
justified heuristically as follows. First of all, the dynamics 
are expected to proceed by ordering not individual grains, 
but entire correlated clusters of typical length £ int . This 
already brings in a factor of £j n t . Furthermore, if we con- 
sider one such piece after a microscopic time but before 
ordering, we see that typical fluctuations of the mean ori- 
entation around its ground-state value (|3.13[) are expected 



We now present quantitative data for the phase di- 
agram sketched in Figure |H The position of the criti- 
cal line £int, c has been obtained by means of the crite- 
rion Kt = 2/3 for large enough systems, as explained in 
the paragraph below (|4.16[) . Figure [7] shows plots of £i n t,c 
against the coupling constant for both shape parameters 
e = 1 and e — 1/<P. In practice, values of £i n t,c larger than 
100 or 200 become very hard to measure with sufficient 
accuracy. 

The critical value £; n t,c is observed to diverge as g — > 1 
from both sides, consistent with the result that mean-field 
dynamics are fast. The fits shown on Figure [7] as dashed 
lines suggest a quadratic divergence of the form 

6„t,c « —^772 (4.20) 

(5-1) 

as mean-field coupling is approached from both sides. The 
amplitude is estimated to be B\ w 11±3, irrespective of e, 
albeit with a relatively large uncertainty. The above for- 
mula will be corroborated below by the finite-size scaling 
law P~2U1) . 

The behaviour of £i n t, c at weak coupling depends on 
whether e is rational or not. For the rational shape pa- 
rameter e = 1, the non-trivial value £i n t,c ~ 28.4 at g = 
|13) is recovered. For the irrational shape parameter 
£ = 1/<P, 6nt, c is observed to diverge as g —> 0. The fit 
shown as a dashed line suggests a linear divergence of the 
form 

U,c « — , (4.21) 
9 

with amplitude Bq rs 0.8. The above divergence can be 
explained in terms of the statics of the model. The qua- 
siperiodic ground state obtained for irregular grains at 
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neither threshold has any effect on the fast (ballistic) or 
slow (activated) nature of the dynamics. 



Fig. 7. Plot of the critical lines for both shape parameters 
e = 1 and e = l/4>. Top: fi n t, c separating phases A and B 
against g for g < 1. Bottom: £i n t, c separating phases A' and 
B against 1/g for g > 1. Dashed lines: fits incorporating the 
divergence laws P~2fJ)) with B% = 11 and (TOTjl with B = 0.83. 



5 = splits into patches of typical length L(g) ~ 1/g at 
small but non-zero coupling. As each of these quasiperio- 
dic patches is retrieved coherently, it is natural to expect 
that the critical interaction length £i n t. c should diverge in 
proportion to the static length L(g) at weak coupling. 

Since the critical interaction length £i n t,c(<7) diverges 
both as g — > and g — > 1 for an irrational shape param- 
eter, it must exhibit a minimum somewhere in the range 
< g < 1. For e = l/'P, this minimal value £i nt , c « 36, 
reached for g ss 0.3, is in the same ball park as the weak- 
coupling value £i n t, c ~ 28.4 for s = 1. Soon after this mini- 
mum is attained, the values of £i n t jC for e = 1 and e = 1/<P 
begin to be nearly identical. One can therefore view the 
minimum in £int.c as a crossover between a phase where 
the dynamics of individual grains (strong shape depen- 
dence) governs the retrieval of weak-coupling-likc ground 
states, and another where the dynamics of clusters (little 
shape dependence) governs the retrieval of mean-field-like 
ground states. This viewpoint makes it clear why the most 
pronounced effects of shape occur before the crossover. 

In the case of the rational shape parameter e = 1, Fig- 
ure [H] shows a comparison between the critical line in the 
9~£int plane and the static and dynamical thresholds g s 
and gd, respectively given by (|3.18[) and (|4.18p . We re- 
call that the dynamics is fully independent of g below the 
dynamical threshold (0 < g < gd), whereas the attrac- 
tors are exactly the dimerised configurations below the 
(larger) static threshold (0 < g < g s ). It turns out that 
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Fig. 8. Plot of notable lines in the £r-£i n t plane for e = 1, 
providing a quantitative version of the left panel of Figure [4] 
Coordinates have been chosen for the sake of clarity. Left full 
line (d): dynamical threshold (|4.18[) . Right full line (s): static 
threshold (I3.18p . Dashed line: mean-field coupling (g = 1). Line 
with symbols: critical line. 



The three lines shown in Figure|5]seem to become close 
to each other in the lower right corner of the plot, i.e., 
at strong coupling, suggesting that one should look more 
closely at this regime. For large values of g, the dominant 
effect is that of the j n field which propagates interactions 
upwards along the column. As a consequence, with our 
choice of fixing the uppermost spin (see (|4.1I) L the system 
will find it more and more difficult to order. This explains 
in qualitative terms why the critical line £i n t,c falls off to 
zero at large g. The data shown in Figure M suggest an 
inverse logarithmic law for £i n t,c> of the form 



* ^ Boo 
Sint,c ~ "; • 

ing 



(4.22) 



For the irrational shape parameter e = 1/^, the data show 
a smooth linear growth with slope l/B^ w 0.275, i.e., 
Boo ~ 3.6. For the rational shape parameter e, the data 
are observed at first to follow those for e = l/'P, and 
then to cross over rather abruptly to a steeper regime of 
growth. The data seem to become asymptotically parallel 
to the static and dynamical thresholds g s and g^. If this 
observation holds quantitatively, we obtain the asymptotic 
slope Boo = 1 for e = 1. 



Vicinity of the mean-field point 

We now turn to the study of the apparently paradoxi- 
cal dependence of the jamming time T on the interaction 
length, in the vicinity of the mean-field point. We have 
encountered another avatar of the same paradox earlier, 
concerning the behaviour (|4.19l) of the velocity V right at 
the mean-field point. For generic values of the coupling 
constant (g ^ 1), the mean jamming time (T) is an in- 
creasing function of £j n t , at least for a large enough system. 
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Fig. 9. Plot of l/£int,c against \ng, for both shape parameters 
e — 1 and e = l/<5. The dashed line has slope 0.275. Upper 
full lines: static (s) and dynamical (d) thresholds for e = 1. 



As shown by (|4. 13f) . f|4.15[) . and fl4.16[) . it grows progres- 
sively as £i n t increases in general: linearly in A" in the bal- 
listic phase, quadratically in A" right at the critical point 
(£mt = Cint.c), and exponentially in A" in the activated 
phase. Right at the mean-field coupling (g = 1) however, 
the divergence law (|4.19l) of the velocity implies that (T) 
is a decreasing function of £j n t . Far from being paradoxical, 
such behaviour is only to be expected; in general (far away 
from the mean-field coupling) , an increase of £j n t implies 
an increase of locally felt frustration, and jamming times 
are thereby increased. Right at the mean-field point how- 
ever, there is in effect no local frustration, simply because 
the grains in the system are effectively non-interacting. 
In its immediate vicinity, an increased correlation length 
£int implies that clusters of interacting grains of typical 
length £j nt reorganise themselves collectively to reach a 
given ground state, so that the jamming time decreases 
with increasing £j nt . 

Figure [TO] shows logarithmic plots of (T) against g/(g+ 
1) for e = 1/<P. The data on the upper panel correspond 
to A~ = 100 and several values of £i n t. The jamming time 
is observed to have a highly non-monotonic dependence 
on g. All the data start from the same value ((T) s» 39) in 
the g — > limit. For £j nt = 1 the data show a broad and 
shallow minimum. As £i n t is increased, the data develop a 
clear minimum near the mean-field coupling (g — 1) and 
a maximum around g « 0.3. The maximum rises suddenly 
for £; nt ~ 35, i.e., near the minimum value £i n t,c ~ 36 of 
the critical line, which is precisely reached for g ss 0.3. 
The observed rise therefore corresponds to the crossover 
between the ballistic and activated phases. The data on 
the lower panel correspond to £; nt = oo, so that the sys- 
tem is in its extreme activated regime, except in the im- 
mediate neighborhood of the mean-field coupling g = 1. 
The observed irregularities are genuine, rather than being 
an artifact due to numerical noise; it was noticed that, if 
the column sizes were chosen to be successive Fibonacci 
numbers, the amount of irregularity could be kept to a 
minimum. 

The data on the lower panel of Figure [TO] show that 
the jamming time exhibits a dip around the mean-field 
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Fig. 10. Logarithmic plots of the mean jamming time (T) 
against g/(g + l) for e = \ j$>. Top: N = 100 and several values 
of £mt- Bottom: £i n t = oo and column sizes equal to Fibonacci 
numbers from F5 = 5 to Fg = 34. Dashed lines: mean-field 
coupling (g = 1). 



coupling (g — 1) in the £j n t = 00 limit, which gets more 
and more symmetric, deep and narrow as A~ is increased. 
The width of this dip can be measured by introducing two 
coupling constants g_ (N) < 1 < g + (N), one on either side 
of the mean-field point, such that (T) = 3 (this value of the 
jamming time is chosen for convenience). Figure [TT] shows 
a plot of the products N 1 / 2 ]ng±(N) against 1/N. Both 
sequences of coupling constants are observed to behave 
very symmetrically, and to depend on N in a very irregular 
way. The hulls of the data however converge to the non- 
trivial limits ±1.75, implying the scaling law 

5g+(N) t* -8g-(N) ~ N~^ 2 . (4.23) 

We have indeed 5g = g — lsslng for Sg <C 1, in conformity 
with the ordinate of Figure [TT] 

The reason for the above scaling behaviour can be ex- 
plained as follows. For £i nt = 00, the local field H n acting 
on grain n reads 

JV 

H n = V -f(a n ) + Sg ]T f(* m ), (4.24) 

m— n+l 

where n is the mean field introduced in (13. 1 1|) . In a random 
configuration of grain orientations, the difference H n — r\ 
can thus be evaluated to be of order (Sg) A -1 / 2 , by the 
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Fig. 11. Plot of the product N 1/2 lng±(N) against 1/N, 
where g±(N) are the two coupling constants such that (T) = 3 
for e = \/$ and £i nt = oo. Upper data: g+(N) > 1. Lower data: 
g-(N) < 1. Dashed lines have intercepts ±1.75 and slopes ±19. 



law of large numbers. This estimate has two consequences. 
First, requiring that the difference H n — 77 is of order unity 
allows one to recover (|4.23l) . Second, the above estimate 
can also be viewed as a finite-size scaling form of the re- 
sult ((47201) . The width 6g ~ N~ 1/2 of the dip is indeed ex- 
pected to be such that £i n t is comparable to N. This yields 
the quadratic divergence £i n t,c ~ 1/ (<5c/) 2 , as in (|4.20p . 

Finally, the irregularities visible in the bottom panel of 
Figure QUI and especially in Figure [TTJ suggest an analogy 
with the phenomenon of defect nucleation, encountered 
in earlier work |12j to explain intermittency in the case 
of a weakly tapped (-T <C 1) column of irregular grains. 
Around the mean- field point, the predominant behaviour 
is of course dominated by the mean field; however, occa- 
sional perturbations in the form of excess j n and h n fields 
can lead to the nucleation of defects at specific sites n, 
generating the self-similar patterns observed. 



Role of a finite £ dyn 

So far we have limited the discussion of the dynamical 
phase diagram to the case where activation energies are 
negligible, so that £dyn = 00. The role of a finite dynamical 
length £dyn has already been investigated in the g — > 
regime for e = 1 [13j . Its main effect is to induce two 
novel dynamical phases, whose main characteristics are as 
follows. 

• Logarithmic phase. The logarithmic phase replaces the 
ballistic one for £; nt < £i n t.e when £dyn is much smaller 
than N. In this regime, the system still orders from above, 
but the growth of the thickness L(t) of the upper ordered 
layer is slowed down by gravity, according to the local 
frequencies (|2.3p . We have therefore 

** F ^hd' (425) 

The results (11712")) and (jiTTBl are recovered for £ dyn > N, 
i.e., in the ballistic phase. In the logarithmic phase, when 



£dyn <C N, the thickness of the ordered layer is predicted 
to grow logarithmically: 

L(t)«£ dyn ln ™ (4.26) 

Sdyn 

so that the jamming time diverges exponentially fast: 

• Glassy phase. The glassy phase replaces the activated 
one when £dyn is sufficiently small. The crossover between 
the activated and glassy phases takes place when the ex- 
ponential growth of the slowest local time scale of the 
column, 1/0jv = e N ^ d y n (which also governs the jamming 
time (I4.27[) in the logarithmic phase), exceeds the entropic 
growth (T) ~ e aN characteristic of the activated phase. 
This line of reasoning predicts that the glassy phase can 
only be observed if £dyn is a microscopic length: £dyn < 
1/a. The jamming time diverges exponentially with N in 
the glassy phase, according to (T) ~ e N ^ dyn . 

Since the dependence on £dyn remains unchanged with 
respect to the earlier case (see [13]), the purely depth- 
dependent features of the glassy and logarithmic phases 
remain the same. Among the novel features of the present 
model, we mention one related to the existence of a mean- 
field point. 

Consider the model right at the mean-field coupling 
(g = 1). In the case of mean-field statics and slow dynam- 
ics, i.e., for £i nt = 00 and N ^> £dyn, the jamming time 
grows exponentially with TV, with an e-dependent prefac- 
tor given by (|4.8p . In the generic situation where the col- 
umn size N is much larger than both lengths £dyn and £i„t , 
the jamming time also grows exponentially with TV, albeit 
with the 'trivial' prefactor of (I4.27[) . Figure H"2l presents a 
logarithmic plot of (T) against the ratio iV/£dyn, for e = 
l/<£, 3 = 1, £dyn = 50, and several values of £i nt . Both lim- 
iting growth rates are clearly observed, the crossover be- 
tween both regimes occurring for very large values of £i n t • 
We note again that the jamming time decreases with in- 
creasing £j„t in the vicinity of the mean-field limit, consis- 
tent with the picture presented earlier. Finally, in order 
to avoid having large irregularities, N has been chosen to 
be a multiple of the Fibonacci number F\q = 55. 



5 Dynamics: attractors 

We next turn to the statistics of attractors sampled by 
zero-temperature dynamics, starting from a random ini- 
tial configuration. This question is interesting because of 
its relationship with Edwards' flatness hypothesis [TB]; ac- 
cording to this, attractors are sampled uniformly, so that 
the static entropy S (see (|3.8[) ) and the dynamical en- 
tropy S (see (|4.10p ) coincide. 

For the time being, we consider the case where £dyn = 
00. Attractor statistics have been investigated in [13] for 
e = 1 in the g — > limit, where ground states have a simple 
characterisation: they are all the dimerised configurations. 
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Fig. 12. Logarithmic plot of the mean jamming time (T) 
against the ratio N/^yu, for e = 1/$, g = 1, £d yn = 50, and 
several values of £i nt . The left dashed line has the unit slope 
corresponding to (|4.27[) . The right one has the slope corre- 
sponding to (H2J), i.e., (1 - e)/(1 + e) = l/i> 3 « 0.236067. 



In the present case, for generic values of the parameters, 
the problem is more difficult because ground states are 
not known a priori, and most certainly do not allow for a 
simple static characterisation. In this respect the present 
situation is similar to that of tapping dynamics on various 
models, e.g. the Kob- Andersen model [21]. In the following 
we present data illustrating what are, according to us, the 
main features of the statistics of attractors. 

First, in order to quantify the role of the coupling con- 
stant g, we have chosen to focus on the dynamical overlap 
between attractors at coupling g and ^ at infinites- 
imal coupling (g = + ), defined as 



Q : 



l 



N 

N ^ 

n=l 



(5.1) 



Note that here the system is started in the same initial 
configuration and subjected to the same stochastic noise 
(i.e., in practice, the same sequence of random numbers). 

Figure [13] shows data for Q for e = 1 and e — l/'P, 
N = 50 and 100, and two values of the interaction length 
(£int = 3 and 10) deep in the ballistic phase, i.e., much 
smaller than £i n t,e- The overlap behaves differently in both 
cases. For the irrational shape parameter e = 1/<P (top), 
the dynamics at infinitesimal coupling drive the system to 
its unique quasiperiodic ground state, so that f2 is noth- 
ing but the overlap between the finite-coupling attractor 

and that ground state. This overlap exhibits a con- 
tinuous decay as a function of g, which is remarkably size- 
independent. This suggests the following picture: quasipe- 
riodic ordering spreads from the top of the column in the 
ballistic regime, as does the splitting of the quasiperiodic 
state into finite patches of length L(g) ~ 1/g at finite cou- 
pling. For a given value of g, the orientation of a grain 
is fixed, depending on its position in one such patch. In- 
creasing the length of the column leaves this orientation 
unchanged, only adding on more, similar patches corre- 
sponding to a given attractor, and thus leaving the over- 



lap function unchanged. For the rational shape parameter 
e = 1 (lower panel of Figure [13]) . the overlap behaves in 
a very different manner. It remains equal to unity in the 
range < g < ga, where the values of the dynamical 
threshold ga (see (|4.18[> ) are shown as arrows. This is ex- 
pected, as the dynamics are strictly independent of g in 
that range. The overlap then falls off very abruptly, more 
and more so for larger systems. 
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Fig. 13. Plot of the dynamical overlap Q against g for both 
shape parameters e — 1/$ (top) and e = 1 (bottom) and 
the same values of fi n t and N. Arrows on top of the lower 
panel: dynamical threshold ga m 0.395612 for £i nt = 3 and 
g d w 0.105170 for £ int = 10. 



We now turn to the statistics of attractors per se. Let 
us first consider in detail the case e = 1. In the weak- 
coupling phase (0 < g < g s ), attractors are known to 
consist of dimers. It is thus natural to characterise attrac- 
tors by their contents in (polarised) dimers all over the 
phase diagram. The main difference between the present 
general situation and that of the weak-coupling regime is 
that now attractors are only partially dimerised. We are 
thus led to introduce two local dimer order parameters: 



Ak = -(C2fc-l(72fc)i 



Ilk = 2^ a2k 



0~2k-l) 



(5.2) 



The first order parameter is such that Ak = 1 whenever 
there is a dimer at the fc-th position and irrespective of 
its polarisation, whereas Ak = — 1 otherwise. The second 
one is sensitive to the polarisation of the dimer at the k-th. 
position; it is such that Ilk = +1 if there is a — f- dimer, 
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Ilk = — 1 if there is a H — dimer, and Ilk — otherwise. 
Both local order parameters vanish on average in a ran- 
dom configuration, where orientations are uncorrelated. 
Assuming that the column consists of an even number of 
grains (N — 2K), we also define global order parameters 
as spatial averages of the local ones: 

fc=l k=l 

The polarisation-sensitive order parameter 77^ can be re- 
cast in terms of the staggered orientation profile 

S n = (-1)'>„) (5.4) 

as 

n k = ^(S 2 k+S 2 k-i), (5.5) 

so that 

1 N 

n = ^s n ( 5 - 6 ) 

n=l 

is nothing but the mean staggered orientation. 

Figure [14] shows a plot of the global dimer order pa- 
rameter A of the attractors against g, for e = 1, N = 20, 
and £j nt = 3 and 10. This modest system size has been cho- 
sen in order to test Edwards' flatness hypothesis fully. The 
plot indeed presents a comparison between (i) the static 
(or a priori) ensemble, where all attractors are obtained 
by means of an exact enumeration and taken with equal 
weights, and (ii) the dynamical ensemble, where attractors 
are sampled according to the dynamics, with a random 
initial configuration. Data pertaining to both ensembles 
behave similarly. They remain equal to unity in the whole 
range < g < g s , where the values of the static thresh- 
old g s (see (|3.18[1 ) are shown as arrows. Both datasets then 
fall off in a similar way (including their fine structure due 
to the finite size of the column), thus indicating that at- 
tractors are sampled rather uniformly by the dynamics. 
In other words, Edwards' hypothesis, although not exact, 
provides a good approximation to the attractor statistics 
in this ballistic regime. 

Next, in order to investigate the possible consequences 
of slow (activated) dynamics on the statistics of attractors, 
we will deal with larger values of £i n t and with arbitrary 
values of e. A natural observable in this case [13] is the 
pseudo-energy £ per grain, defined as 

1 N 

£ = --Y, H n°n- (5.7) 

n=l 

This definition can be motivated as follows. If the a n were 
independent spins in external fields H n , (|5.7p would be 
the corresponding Hamiltonian. In the present model, the 
local fields H n depend on the orientations a m in a non- 
symmetric way, so that the dynamics do not obey detailed 
balance, and the statics are not described by a Hamilto- 
nian. The pseudo-energy defined by (|5.7|) however provides 
a useful measure of the amount of disorder, either in an 
arbitrary configuration or in an attractor. 
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Fig. 14. Plot of the global dimer order parameter A of attrac- 
tors against g for e = 1, N = 20, and two values of £i n t- Lines 
with symbols: dynamical ensemble. Lines without symbols: 
static ensemble. Arrows on top: static threshold g s ~ 0.692394 
for Cint = 3 and g s « 0.473376 for £ int = 10. 



We have chosen to focus on two values of the coupling 
constant, one on each side of the mean-field coupling, i.e., 
g = 0.6 and g = 2, where the critical value of the inter- 
action length is respectively large and small. These are 
£int,c ~ 51 for £ = 1 and 48 for e = 1/$ when g = 0.6; 
and £mt, c ~ 5.2 for e = 1 and 5.9 for e = 1/& when 
.9 = 2. Figure [TBI shows plots of (minus) the mean pseudo- 
energy {£) per grain against £i nt in a range containing the 
critical values £i n t,c (shown by arrows) which separate the 
ballistic and the activated phases. The mean energy ex- 
hibits a very weak and regular increase as a function of the 
interaction length £i nt ; this verifies what we might expect, 
that increasing correlations will increase the amount of 
order in the system. The mean energy however, shows no 
anomalies at all as the critical point is crossed; this smooth 
behaviour is to be contrasted with the explosive rise in 
jamming times which accompanies it. Over the ranges of 
values of £i n t corresponding to the plotted data, the mean 
jamming time indeed increases by factors of 150 for e = 1 
and 580 for e — 1/<P in the case of g = 0.6, and by factors 
of 2100 for e = 1 and 80 for e = 1/<P in the case of g = 2. 

Role of a finite £d yn 

The effect of a finite dynamical length £dyn on the statis- 
tics of attractors has already been investigated in [13] for 
e = 1 and g —> 0, that is, for weak frustration. The main 
qualitative conclusion there was that a truly non-trivial 
sampling of attractors was observed only in the glassy 
phase, whereas the attractor statistics in the three other 
dynamical phases were found to be in qualitative agree- 
ment with Edwards' flatness hypothesis. The effect of in- 
creasing frustration, in the present paper, might be ex- 
pected only to enhance the non-triviality in the sampling 
of attractors in the glassy phase, thus deepening the con- 
trast between this and the other three phases. 

We consider first the vicinity of the mean-field point. 
Figure [L7J1 shows the orientation profile (a n ) of the attrac- 
tors for parameters similar to those used in Figure 1121 
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Fig. 15. Plot of (minus) the mean attractor pseudo-energy {£) 
per grain against £i nt , for N = 200 and both shape parameters. 
Arrows: critical values £i n t,c- Top: g = 0.6. Bottom: g — 2. 








Fig. 16. Plot of the orientation profile of the attractors again 
depth n for N = 100, e = 1/$, at the mean-field coupling 
g = 1, and £dyn = 50. Symbols: data for several values of £i n t- 
Dashed line: mean orientation (a) = — 1/$ 3 ~ —0.236067. 



i.e., e = l/#, g = 1, N = 100, and £dyn = 50. In the case 
of mean-field statics (£i„t = oo), a non-trivial orientation 
profile is observed, similar to that shown in Figure [3] In 
the generic situation where the interaction length £; nt is 
finite, the profile soon becomes very nearly flat, and there- 
fore equal to its mean value (|3.5[) , i.e., (a) = — 1/4> 3 ss 
—0.236067, thus indicating that attractors are sampled 
nearly uniformly by the dynamics. 

Next, we explore the statistics of attractors for a gene- 
ric coupling g 7 to, in particular, test the validity of Ed- 



wards' hypothesis. We use a coupling constant greater 
than the mean field value (g > 1), so as to be able to have 
a relatively small £i n t,c with our choice of boundary condi- 
tion l|4.1[) . The consequence of this is that the slow phases 
(activated and glassy) are easily observed with a modest 
system size. This possibility of tuning £i n t,c, and conse- 
quently system sizes, to reasonable values was not present 
in the case of the weak-coupling regime [13] , where lengths 
were perforce large; it explains in part our rationale for not 
changing the boundary condition (|4.1[) , of which more will 
be said in the concluding section. 

We choose e = 1, in order to keep using the local dimer 
order parameters (|5.2p . a coupling constant of g = 2, so 
that £int, c ~ 5.2, and a system size N — 50. We compare 
the nature of the attractors at £j n t = 2 (below £i n t. c : bal- 
listic to logarithmic crossover) and £i n t = 8 (above £i n t, c : 
activated to glassy crossover) . Figures [17] to [19] respec- 
tively show plots of the staggered orientation profile S n 
and of the local dimer order parameter Ak and 77^ for 
several values of £dyn spanning the crossovers. 
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Fig. 17. Plot of the staggered orientation profile of the at- 
tractors against depth n for e = 1, g = 2, N = 50, and several 
values of £d yn - Top: £i nt = 2. Bottom: £i nt = 8. 



Our earlier speculations are concretised by the follow- 
ing observations. The data for £; nt = 2, shown in the upper 
panels of Figures [T7] to [T51 and pertaining to the crossover 
between the ballistic and the logarithmic phases, hardly 
show any dependence on £dyn, i-e., on the fast (ballis- 
tic) or slow (logarithmic) nature of the dynamics. Since 
we have established that Edwards' flatness holds, at least 
qualitatively, for the ballistic phase, this is a clear indica- 
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Fig. 18. Same as Figure [T71 for the dimer order parameter Ak, 
plotted against the dimer number k. 
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Fig. 19. Same as Figure[T7Jfor the polarisation-sensitive dimer 
order parameter Ilk, plotted against the dimer number fc. 



tion that it holds also for the logarithmic phase. In terms 
of the actual ordering, we note the high and nearly uni- 
form value of the dimer order parameter (Ak ~ 0.75) in 



the upper panel of Figure [T5J the picture is that the at- 
tractor is most ordered at the base, and more or less or- 
dered (dimerised) throughout its length, except near the 
very top. The staggered orientation profile and the order 
parameter Ilk, shown in the upper panels of Figures [171 
and [TDJ demonstrate that the dimers are essentially un- 
polarised, except in a thin boundary layer near the top of 
the column. 

On the other hand, the data for £j nt = 8, shown in the 
lower panels of Figures [TTJ to \T9\ and pertaining to the 
crossover between the activated and the glassy phases, 
exhibit a strong dependence on £dyn- Highly non-trivial 
profiles are observed for the smaller values of £dyn, and 
especially for £dyn = 5, which can be considered as glassy, 
given the modest size of the column. In this case, we see 
that the ensuing structure of the column is highly hetero- 
geneous. The large values of the parameters Ak and Ilk 
observed in the upper part of the column indicate a large 
degree of dimerisation, most of the dimers being polarised 
as — h In other words, there is a definite preference for 
one of the 'crystalline' arrangements of dimers found in 
earlier work [T3] in the glassy regime for e = 1 and g -C 1, 
The choice of one preferred direction of polarisation can 
be explained as follows. The boundary condition <j\ = +1, 
together with the observed smallness of (02), yields a def- 
initely negative /13, and hence a trend toward (0-3) < 0, 
i.e., II 2 > 0. For £ int = 8 and £ dyn = 5, (a 2 ) « 0.04 is 
indeed much smaller than (173) ~ —0.62. The same po- 
larising effect acts on the deeper dimers as well. Another 
notable difference with respect to the weak-coupling situa- 
tion considered in [T3] is that the lower part of the column 
appears much less ordered. One may speculate that a dis- 
ordered lower part is a generic characteristic of attractors 
in the glassy regime. Grains are indeed very slow as soon as 
n £dyn, and thus hardly equilibrate. The weak-coupling 
regime for e = 1 appears as an exception to this general 
rule, because there it is already known from statics that 
all ground states are fully dimerised. 



6 Discussion 

The full phase diagram of the frustrated column model 
has been presented in this work. This model of a column 
of grains has been developed [TDJ and investigated in ear- 
lier work, first in the directed situation (g — 0) for arbi- 
trary values of the shape parameter e [11112] , and then in 
the weak-coupling regime of the symmetric case (e = 1, 
g <C 1) [13]. The present work is the first investigation of 
the model over its entire parameter space, i.e., for generic 
values of e, g, and of both lengths £i nt and £d yn - One of the 
most novel features is the existence of a mean-field point 
(g = 1) where the erstwhile local constraints on grain ori- 
entations disappear, becoming global; the model in this 
limit is similar to one of non-interacting grains, which has 
been analysed in [ID] . 

The case of e = 1 provides a useful illustration of many 
features of the model, including the physical nature of the 
mean-field limit. In this symmetric situation there are es- 
sentially two ways in which the average orientation of zero 
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can be achieved: by dimerised packings in the presence of 
strong local compacting constraints (g — > 0), and by uni- 
formly random packings such that there are equal numbers 
of ordered and disordered grains globally, in the absence of 
local constraints (g — 1). Most of the features of the phase 
diagram for e = 1 for g < 1 can be explained in terms of 
an interpolation between these two extremes. Beyond the 
mean-field limit (g > 1), the picture is one of a column of 
grains that is strongly frustrated, partly due to the cho- 
sen boundary condition (the uppermost spin is fixed), so 
that the activated phase extends down to arbitrary small 
values of £i„t at strong coupling. 

This is a natural point at which to discuss the implica- 
tions of our choice (|4.ip of boundary condition. For the di- 
rected model [11112] . where order propagated downwards, 
it was a natural choice to fix the top spin. This choice 
may seem to become less and less natural as increasing g, 
as the effects of the frustrating field j n made the propa- 
gation of order less and less directional. Indeed, it might 
be argued that, in the large- coupling regime, the natu- 
ral choice would be to fix the bottom spin, corresponding 
to the predominant upward propagation of order. There 
are two aspects to this issue. On the one hand, fixing the 
bottom spin rather than the top one might in fact make 
the phase diagram look more symmetric with respect to 
the weak-coupling end, and in particular show the explicit 
reappearance of the difference between regular and irregu- 
lar grains in this. On the other hand, this would imply that 
we were replacing a fully frustrated column with a column 
where, instead, the 'reverse' field j n could yield ballistic 
propagation until a finite £i n t,c w &s reached - which is not 
physically correct. Thus our choice of fixing the top spin 
throughout the phase diagram has the advantage of re- 
taining the sense of the interpolation from a fully directed 
model without frustration to the fully frustrated model 
for large g. 

What can we predict experimentally for a box of gra- 
ins? First, the ubiquity of our four phases (ballistic, loga- 
rithmic, activated and glassy) throughout the phase dia- 
gram (except at the mean-field point) vindicate our earlier 
picture [S] that the top of such a box would look ballistic, 
the middle activated and the bottom glassy. The relative 
sizes of the phases would of course vary depending on 
frustration: we might expect, given our study, that in the 
presence of strong frustration the size of the glassy phase 
would increase, and that of the ballistic one decrease. For 
timcscales that are typical of experiment or simulation, we 
would see the highest fluctuations coming from the acti- 
vated phase in the middle of a typical box [S]; however, for 
much longer times of observation, we would see large non- 
ergodic fluctuations at the bottom of the box, consistent 
with the glassy phase. 

Finally, we predict that the effects of shape are most 
likely to be visible for very weak or very strong frustra- 
tion, where there is a dominant propagation of order in 
our column in a given direction. In between them, when 
there is a balance between the propagation of order due 
to gravitational settling and that due to frustration, we 



might expect shape effects to disappear as clusters of gra- 
ins became the units of reorganisation. 
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